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Abstract 

New schemes for the nuclear spin quantum memory are proposed based on a system composed 
of two electrons or one electron coupled to a single nuclear spin in isotopically purified group IV 
elemental and II- VI compound semiconductors. The qubit consists of the singlet state and one of 
the triplet states of two electrons or simply of an electron spin. These systems are free from the 
decoherence due to the nuclear dipole-dipole interaction and are advantageous for the long memory 
time. In the case of two electrons, the protocol for the quantum state transfer between the electron 
spin qubit and the nuclear spin qubit is based on the magnetic or electric field tuning of the singlet- 
triplet state crossing and on the hyperfine coupling supplemented with a well-defined scheme to 
initialize the nuclear spin. In the case of a single electron qubit, the quantum state transfer is 
driven by the hyperfine interaction itself without the need of the nuclear spin initialization. Many 
practical systems are considered, e.g., two electrons loaded on a Si or ZnSe quantum dot, a single 
electron charged state in a Si quantum dot doped with a P atom, a single electron charged 28 Si 
quantum dot doped with an isotope atom of 29 Si, and a localized electron system of Si:P and 
ZnSe:F in the bulk crystal. General aspects of these systems are investigated and a comparison of 
merits and demerits is made between the two-electron qubit and the single-electron qubit. 
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I. INTRODUCTION 



Recently, a quantum media converter from a photon qubit to an electron spin qubit was 
proposed for quantum repeaters^. Quantum information can take several different forms 
and it is preferred to be able to convert the information among different forms. One form 
is the photon polarization and another is the electron spin polarization. Photons are the 
most convenient medium for sharing the quantum information between distant locations. 
However, it is necessary to realize a quantum repeater in order to send the information 
securely over a very long distance overcoming the photon loss. A quantum repeater requires 
two essential ingredients, namely, the quantum state transfer between a photon and an 
electron spin^ and the correlation(Bell) measurement between two electrons created by the 
quantum state transfer from two different photons^. Furthermore, it is desirable to have a 
quantum memory to store the quantum state of the electron spin. The most suitable medium 
for that purpose is the nuclear spins because of their extremely long coherence time. It has 
been already proposed to use the collective nuclear spin state of the host medium, e.g., GaAs, 
as a quantum memory, employing the hyperfme (hf ) interaction to transfer the electron spin 
state to the nuclear spins^. However, in that proposal one has to achieve a high degree of 
nuclear polarization for the quantum memory of high fidelity, which has been prohibitive so 
far due to the low nuclear spin polarization achievable experimentally^. In III-V compound 
semiconductors, all the nuclei have a nonzero spin and the nuclear spin quantum memory 
is necessarily subject to the decoherence induced by the nuclear dipole-dipole interactions. 
Thus it is advantageous to employ a system with a few nuclear spins like the nitrogen- 
vacancy (NV) center in diamond for the quantum memory or register. The nitrogen spins 
of NV centers in diamond^ and phosphor spins in Si:P— have been demonstrated to be 
promising as quantum memories. In these systems the degeneracies associated with nuclear 
spins are lifted by the hf interaction, enabling a selective addressing of nuclear spin states 
by external microwave or radio frequency (rf) fields. 

In view of these progress, we consider isotopically purified materials made of the group 
IV elemental or II- VI compound semiconductors, such that the number of atoms with a 
nonzero nuclear spin can be reduced down to only one and propose a few new schemes for 
the nuclear spin quantum memory in which the hf interaction itself drives the quantum 
state transfer (QST) between the electron spin qubit and the nuclear spin qubit. We study 
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a system composed of two electrons or one electron and only one nuclear spin. Actual 
examples for the two-electron qubit are two electrons loaded in a 28 Si quantum dot (QD) 
doped with an isotope atom of 29 Si or in a ZnSe QD doped with an isotope atom of 77 Se and 
a single-electron charged state in a Si QD doped with a P atom or in a ZnSe QD doped with 
a fluorine (F) atom. In these examples, the qubit is composed of the singlet state and one 
of the triplet states of two electrons and the QST is operated at the singlet-triplet crossing 
point. This crossing point can be approached by tuning the magnetic field or the electric 
field. The QST between the electron spin qubit and the nuclear spin qubit can be carried 
out by the hf interaction itself, reinforced with a well-defined scheme for the nuclear spin 
initialization which is based on the electron spin state measurement^. A key requirement 
in this QST is that the singlet-triplet anticrossing gap should be much smaller than the hf 
interaction energy. We reveal that this requirement is satisfied favorably in the donor-bound 
electron system but not in the delocalized electron system in a QD. 

On the other hand, for the single-electron qubit, practical examples are a single-electron 
charged 28 Si QD doped with an isotope atom of 29 Si, a single-electron charged ZnSe QD 
doped with an isotope atom of 77 Se, and a localized electron system of Si:P and ZnSe:F in 
the bulk crystal. Here the electron spin plays the role as a qubit and the single nuclear spin 
of an isotope atom or a donor atom plays the role as a quantum memory. It is advantageous 
that the relative magnitude between the spin-orbit interactions and the hf coupling energy 
is not relevant in this QST. A more important feature in this QST is that the nuclear spin 
initialization is not necessary. 

The paper is organized as follows. We first present the results on the magnetic field 
tuning of the singlet-triplet state crossing of two electrons loaded on a single QD and of two 
electrons in a single-electron charged QD doped with a donor atom. We discuss the effects 
of the spin-orbit interactions which induce the singlet-triplet state anticrossing and become 
obstacles to our QST protocols. Then we discuss the hyperfine interaction and propose a few 
protocols for the QST between the electron spin qubit and the nuclear spin qubit. We also 
investigate the system composed of a single electron and a single nuclear spin and present 
a QST protocol between them, because the single-electron qubit is more fundamental as 
a building block of devices for the quantum information processing. Finally, we make a 
comparison of merits and demerits between the two-electron qubit and the single-electron 
qubit and conclude that the donor-bound single electron system is preferable with respect 
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to the fast QST operation, no need of the nuclear spin initialization and the irrelevance to 
the relative magnitude between the spin-orbit interaction and the hyperfme interaction. 



II. MAGNETIC FIELD TUNING OF THE SINGLET-TRIPLET STATE CROSS- 
ING OF TWO ELECTRONS IN A QUANTUM DOT 

A. Energy spectrum of a pair of electrons in a Si or ZnSe quantum dot 

First we consider a pair of electrons in a Si QD. In bulk Si, the conduction band minima 
have a six- fold degeneracy, which makes it difficult to define a robust electronic qubit. 
However, this degeneracy is lifted in Si quantum well (QW) structures. The strain in the 
lateral direction (x and y) lifts up the four conduction band minima along the x and y 
directions by about 100 meV— . Additionally, the confinement in the growth direction couples 
the lowest two conduction band minima along the z direction, removing the degeneracy 
completely^. At the conduction band minima along the z direction, the Bloch functions, 
which are separated by an energy difference about 1.5meV— , are given by 



where (0, 0, ±/co) {ko — 0.857r/a) is the wave vector at the band minima along the z direction 
with the lattice constant a of Si, <fi a phase factor related to the valley mixing, and Wfc (r) 
is the periodic part of the Bloch function at the band minima normalized in the unit cell 
volume. The lower-energy one of ip+ and ip_ is determined depending on the QW thickness^. 
The actual wavefunction is given by the product of the lowest energy Bloch function and 
the envelope function F(r) satisfying the effective mass equation. Under a magnetic field 
along the z direction the effective mass equation for the envelope function is given by 



where U represents the harmonic confinement in the lateral direction, V is an additional 
potential in the z direction describing, e.g., the electrostatic confinement, m t (= 0.19 mo,mo 
being the free electron mass) and mi(= 0.92 mo) are the transverse and longitudinal effective 



^+(r) = y/2cos(k z + <f))\u ko (r)\ , 
^-(r) = V2sm(k z + (f>)\u ko (r)\ , 



(1) 
(2) 
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masses, respectively. The Zeeman energy term does not appear in Eq. ([3]) because F(r) 
represents only the orbital part. But the Zeeman energy is included in the calculation of 
the energy level structure. 

For a vanishing magnetic field, the ground state of two electrons is the spin singlet stated. 
However, in the presence of a magnetic field, the energy spectrum does not necessarily 
follow the Lieb-Mattis theorem. The magnetic confinement and the Coulomb interaction 
together lead to crossings in the energy spectrum. In FigJH the exchange energy J = 
^(triplet) — .E(singlet), namely, the energy difference between the ground singlet state and 
the excited triplet state of two electrons, is plotted for two values of the lateral confinement 
energy fkjQ. The point at J = indicates the singlet-triplet crossing point. When electrons 
are confined weakly in the lateral direction, the effect of the magnetic field on the orbital 
motion begins to appear at the weak field, influencing the exchange energy between two 
electrons. Conversely, under the strong confinement in the lateral direction, the effect of the 
magnetic field on the exchange energy between two electrons becomes manifest at the strong 
field, shifting the singlet-triplet crossing point to a higher magnetic field. These qualitative 
features can be confirmed in Fig. 1. 

Now we discuss the effect of the additional potential in the z direction produced by, e.g., 
electrical gates. To facilitate the arguments, that effect is taken into account by assuming a 
harmonic potential for V(z) in Eq. ([3]) as 

V(z) = ~m z u, 2 V , (5) 

where hu z is the harmonic confinement energy. In general, the spatial extension in the 
growth direction enhances the three-dimensional character of the electron motion and leads 
to the reduction in the Coulomb energy and to the weaker dependence of the orbital motion 
on the magnetic field, pushing the singlet-triplet crossing point to higher magnetic fields. 
In FigfSJ the effect of the additional potential on the energy spectrum is exhibited. As 
the confinement energy fiw z in the z direction is increased, the singlet-triplet crossing point 
moves to higher energies because of the increase in the Coulomb energy and shifts to weaker 
magnetic fields because the two-dimensional character of the electron motion is enhanced 
and the orbital motion becomes more susceptible to the magnetic confinement. 

As an example of the direct-gap material, we consider a ZnSe QD. Contrary to Si, the 
bulk ZnSe has the conduction band minimum at the T point with an isotropic effective mass. 
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FIG. 1: The exchange energy J = ^(triplet) — (singlet), namely, the energy difference between 
the ground singlet state and the excited triplet state in a circularly symmetric quantum dot is 
plotted as a function of the magnetic field strength. The intercept with the line J = represents 
the singlet-triplet crossing point. The solid (dashed) line corresponds to the Si (ZnSe) quantum 
dot. Curves on the left (right) are plotted for hw$ = 1(5) meV. 

The conduction band electron can be described by the wavefunction: 



where Uq(t) is the Bloch function at the T point and the envelope function F(r) satisfies Eq. 
(J2J) with a modification of m t = mi. Thus the physics of the singlet-triplet crossing is the 
same as discussed above. The crossing behaviors are shown in FigfTJ employing the effective 
electron mass m t = mi = 0.16 mo and the dielectric constant k = 9.1—. 

The singlet-triplet crossing of a pair of electrons on a GaAs QD was observed by the gate 
voltage tuning^. The gate voltage controls the number of electrons in the QD as well as 
the shape of the QD potential electrostatically, enabling the ground state tuning. In Fig. [3l 
the dependence of the exchange energy J on the lateral confinement energy ftwo is depicted 
for a few values of the magnetic field strength. When the lateral confinement is increased, 
the orbital energies increase directly proportional to the confinement energy. At the same 



tt(r) = u„(r) F(t) 



(6) 
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FIG. 2: Energy level crossing between the ground singlet (growing curve with the magnetic field) 
and the triplet (decreasing curve with the magnetic field) states in a Si QD with the fixed lateral 
confinement energy fkuo = 1 meV and the additional confinement energy in the growth direction, 
huj z = (a) 1 (solid curves), (b) 5 (dashed curves), (c) 25 meV (dot-dashed curves) and (d) oo 
(dotted curves). As ftw z is increased, the level crossing occurs at a higher energy and at a smaller 
magnetic field. In this plotting the confinement energy in the z direction is subtracted for the ease 
of comparison. 

time, the spatial overlap between electron orbitals is enhanced and consequently the direct 
and exchange Coulomb energies increase. This leads to the increase of J in proportion to 
the lateral confinement energy huo, as seen in Fig. [31 The lateral confinement energy or the 
potential curvature can be controlled by tuning the voltages on several gates. Thus, either 
by electrically modifying the lateral confinement potential or by changing the magnetic field, 
one can tune the singlet-triplet crossing. 

B. Effects of spin-orbit interactions 

So far we have neglected the effect of the spin-orbit (SO) coupling. But this effect should 
be examined because the SO coupling mixes the spin states through the orbital state mixing 
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FIG. 3: The exchange energy J, i.e., the singlet-triplet energy difference, is plotted as a function 
of the lateral confinement energy ftujQ for the magnetic field strength B = 0.4, 0.7 and 1 Tesla from 
top to bottom. 

and affects the singlet-triplet crossing^. If the spin-orbit interaction has a non-zero matrix 
element between the singlet (S) and triplet (T) states, the level crossing turns into the level 
anticrossing which becomes an obstacle to the QST driven by the hf interaction. In order 
to achieve securely the electron-nuclear QST, it is required that the matrix element of the 
spin-orbit interaction is much smaller than that of the hf interaction. In the following we 
will examine this requirement. 

The spin-orbit coupling for the conduction band electron in the linear approximation 
with respect to the momentum operator is given as 

Vso = a R (a x p y - <j y p x ) + a D (a x p x - a y p y ) , (7) 

where the first term is the Rashba SO coupling due to the structural inversion asymmetry^ 9 - 
and the second term is the Dresselhaus term arising from the bulk inversion asymmetry^. 
Here the momentum p is meant by the gauge invariant kinetic momentum —ihV + eA/c 
with the vector potential A, if the system is under a magnetic field. In order to discuss 
the S-T anticrossing, it is convenient to eliminate the SO coupling by applying a unitary 
transformation^- to the Hamiltonian as discussed in Appendix A. The original single electron 
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-P 2 + a R (a x p y - a y p x ) + a D (a x p x - a y p y ) + -g e \i B B 



(8) 



Hamiltonian is given by 

where g e is the electron g-factor. By rotating the coordinate system by an angle tt/4 in the xy 

plane, namely, by introducing a new coordinate system (£, 77, z) defined by = (1, 1, 0) / v2, 

e,j = (—1, 1, 0)/y/2, e z = (0, 0, 1), the Hamiltonian is rewritten as 

1 I 
H = 2"^P 2 - P£cr v (a D + a R ) + p v a^(a R - a D ) + -g e /i B (B^ + ^ct^ + B z a z ) 



2m 



Pi 



h 



-cr, 



2m 



2m 



+ ^9e^Bi.B^a^ + S^cr,, + B z cr z ) - m(a|, + a|) 



with 



A g A 



h 



(9) 

m(a R + a D ) ' " v m(a R - a D ) ' ^ 
where the last constant term in Eq. ([9]) will be omitted. Now, in order to eliminate the 
original spin-orbit coupling, we apply the unitary transformation: 



H=U ] HU with U = exp^ov, - i^-tr^ . 



At 



A r 



(11) 



Then we have 
H 



2m 

h 



P 2 + i^9e^B{B^a^ + B v cr v + B z a z 
— a € + ?-<tA B 7 - ( $-B f + ■?-£ 



A 



At 



Ar? ^ 



(7- 



-L z a z + ^(-2£ 2 p„ + T^pJ) + ^(2?7 2 ^ - £{77,pJ) 
A 5 ™n 



;i2) 



where {A, 5} = AB + Z?A, the first line represents the single electron Hamiltonian under a 
magnetic field without the SO coupling, the second line the SO induced Zeeman interaction 
H§° and the third line stands for the renormalized SO interaction H ren . The energy level 
structure of H is exactly the same as that of the original Hamiltonian H. Thus we can 
discuss the S-T crossing behavior based on the transformed Hamiltonian H. After this 
unitary transformation, the original SO interaction is eliminated but new terms appear with 
a smallness parameter defined by e = it/Xso which is typically about 10~ 3 , where It is the 
lateral extent of the electron wavefunction and A50 represents symbolically the SO length 
A^ and A^. Among the newly appeared terms, the SO induced Zeeman interaction H§ is 
the most dominant term which is of the first order in the smallness parameter e. However, 
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the matrix element of this term between the singlet state IS*) and the triplet states \T±), i.e., 
(S\H§°\T±) can be made to vanish by tuning the direction of the magnetic field as shown 
in Appendix A. When one of the SO coupling constants ap and a# is much larger than the 
other, the perpendicular (z-directed) magnetic field is most favorable. However, the actual 
magnitude of the S-T anticrossing gap is determined by the higher order perturbation terms 
of H§ and contributions from the renormalized SO term H ren . The typical magnitude is 
about 1CT 12 ~ 1CT 11 eV as shown in Appendix A. 

So far we discussed the linear-in-momentum SO coupling. However, the cubic-in- 
momentum SO term is present in general and more detailed arguments are necessary. As 
discussed in Appendix A, the Dresselhaus SO term can be decomposed as 

V D = y£> + V® (13) 
with = {<j>{z)\V D \<j>{z)) = 1 ((]){z)\k' i z \cj){z)){-k x G x + k y a y ) , (14) 

v d ] = l(k x k 2 y a x - kykla y ) , (15) 

where <p(z) is the ground state orbital in the ^-direction, the linear term is already 

(3 s ) 

included in the above Vso an d denotes the cubic-in-momentum SO term. Unfortunately, 
this cannot be eliminated by the unitary transformation as shown below 



~2i 2i 

—(-h{pt, Pr) } + z({p£,p„} - 2p 3 v )0 + T-(-^{p S ,pJ + i{{p%Pn) ~ 2p|)77) 
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<y z + 



Ah 3 

= y ( D + y§ )ren + ■•• , (is) 

where the newly appeared terms denoted by V^ ren are smaller in magnitude than the 
original terms due to the smallness parameter e = it/^so ~ 10 -3 . Furthermore, V^ ren 
does not contribute to the matrix element between the singlet state \S) and the triplet 
states \T±), because this does not change the magnetic quantum number. Thus we have to 
discuss the effect of on the S-T anticrossing gap . As shown in Appendix A, the matrix 
element (S\V^ \T±) vanishes, if the orbital excited state associated with the T± is chosen as 

e rd( x > y) = yi. x T iy) 9rei{x, y) with g rei {x, y) = —= exp[--^ (x 2 + y 2 )} . (17) 
it tty/TT li t 

The higher order perturbation terms concerning V^ 3 "* and the contributions from Vp Ten 
are estimated in Appendix A. The typical magnitude of the S-T anticrossing gap is about 
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10~ 12 ~ 10~ n eV, which is of the same order of magnitude as that for the linear-in- 
momentum SO coupling. Consequently, we can summarize that the S-T anticrossing gap 
is comparable to the hyperfine coupling energy as discussed in Sec. IV for the case of two 
delocalized electrons in a QD with an isotope atom of the host material. This situation is 
unfavorable for our protocol of the nuclear spin quantum memory which requires that the 
S-T anticrossing gap is much smaller than the hyperfine coupling energy. 

However, we consider also the case of the single electron charged QD with a donor impu- 
rity having the nuclear spin in the next Section and find that the S-T anticrossing gap can 
be much smaller than the hyperfine coupling energy. Thus our protocol of the nuclear spin 
quantum memory will be fully effective in this case. 



III. ENERGY SPECTRUM OF A SINGLE ELECTRON CHARGED QD DOPED 
WITH A DONOR ATOM 

Here we investigate the energy spectrum of a singly charged QD doped with a donor. 
The extra electron will hybridize with the donor electron so as to form the singlet-triplet 
states. We can envision three regimes as shown in Fig. HI namely, (a) two electrons are 
bound to the donor, (b) one electron is bound to the donor and the other is delocalized 
over the QD, and (c) both electrons are delocalized over the QD. These regimes depend on 
the position of the donor atom. When the donor is doped at the central part of the QD, 
two electrons are bound strongly by the donor atom, whereas when the donor is doped at 
the peripheral region of the QD, two electrons are trapped weakly by the QD confinement 
potential. In the intermediate case, the donor atom and the QD potential provide the same 
order of confinement for a trapped electron and there occurs the situation that one electron 
is bound to the donor and the other is delocalized over the QD. For cases (a) and (c), it is 
already known that there occurs a singlet-triplet level crossing induced by the magnetic field 
and the Coulomb interaction.—^^ In this section, we focus on the case (b) and confirm that 
the singlet-triplet level crossing occurs at appropriate values of the distance of the donor 
atom from the center of the QD. 

The Hamiltonian for two electrons in a negatively charged donor-doped QD is given by 

H = fe(n) + h(v 2 ) + e2 . (18) 

k\ti — r 2 | 
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(a) 




(c) 



FIG. 4: Schematic configuration of two electrons in a single electron charged QD doped with 
a donor atom. The large (small) parabola represents the QD confinement potential (Coulomb 
potential by an ionized donor). A black dot indicates an electron. Three regimes are depicted: 
(a) two electrons are bound to the donor, (b) one electron is bound to the donor and the other is 
delocalized over the QD, and (c) both electrons are delocalized over the QD. 



Here h is the single-particle Hamiltonian given by 



h -- 

with T 

u 



T + hz + V a + V b 
1 , e . , 



2m 



c 



^g e fi B Ba z , V a = -m(u%p 2 + uj 2 z z 2 ) , 



K\r 



dx\ 



P 



x 



y 



A 



B 



-y,x,0) 



(19) 



(20) 



where T is the kinetic energy part including the effect of a magnetic field applied along the 
z-direction, hz the Zeeman energy with the electron g-factor g e and the Bohr magneton 
/is, V a the confinement potential by the QD, V& the Coulomb potential by an ionized donor 
displaced by a distance d from the center of the QD and x is the unit vector in the x-direction. 

In the Heitler-London approach, the orbital wavefunctions for the spin-singlet (^+) state 
and the spin-triplet (^_) state are approximated as 

1 



with S = 



(0 a (ri)0 6 (r 2 ) ± a (r 2 )0 6 (ri)) 



(21) 
(22) 
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where a and <\> b are appropriate single-particle wavefunctions specified later. Then we have 

E± = ($?±\H\y±) (23) 
(e a + e b + (4> a \V b \(f> a ) + (4> b \V a \(f> b ) + V d ) 



1±\S\ 



± Y^r^ii^lhl^S* + (falhWJS + V x ) 



12 

,2 



with V d = (0 a ( ri )0 6 (r 2 )|^ r |0a(ri)0 6 (r 2 )) , (25) 

ac|Ti — r 2 | 

V* = <^a(ri)0 6 (r 2 )|^-^ , |0 a (r 2 )0 fe ( ri )) , (26) 

ft| r l — r 2| 

e a = (<f> a \T + V a \<f> a ) , e b = (<f> b \T + H|^> , (27) 

where V^(Vx) is the direct (exchange) Coulomb energy. When only the <j) a can be chosen to 
be an exact eigenstate: 

(T + V a )<f> a = e a a , (28) 

the exchange energy, namely, the energy difference between the singlet and the triplet states, 
is calculated as 

21 9I 2 

J = E- — E + = - e ) + (0 a |Wa) + + ^) 

- 1 _ 2 |g|4 ({0alHl0 6 )^ + ( W|0a)S + Fx) + ge^m, , (29) 

where the Zeeman energy is included with the magnetic quantum number m z of the triplet 
state 

Now, (p a is chosen as the exact ground state in the QD confinement potential V a : 

11 o 2 z 2 

(30) 



with ^=V^' / *=v^' n= v^ + T' w -= 1 ^' (31) 

where cj c is the cyclotron frequency and the eigenenergy e a is + hw z /2. On the other 
hand, the donor-bound eigenstates are calculated variationally in the Gaussian functions to 
facilitate the following calculations. The ground state is approximated as 



^ (r) = ^^| eXph 24-24 ] ' (32) 
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where a and (3 are the variational parameters and the corresponding energy is calculated as 

„2 



E\ 



IT- 



k r 



" ~ 2 ~ a 



In 



for a > (3 , 



with R 



me 



2! _ 4 /s\ for /5 = a , 

-^01= arctan [ - 1] | for (3 > a 

H 2 K f\jjJr 



o as 
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(33) 



(34) 



(35) 



H 2 k 2 " me 2 1 1 2R 
In the absence of a magnetic field there is no preferential direction in the donor-bound 
state and thus a = (3. Under a magnetic field along the ^-direction, the squeezing of the 
wavefunction occurs in the transverse direction and a> (3. 

The eigenstate of an electron in the presence of a magnetic field and an ionized donor at 
a position shifted from the origin is determined by 



T 



(36) 



n\r — dx\ . 

The shift of the position of the donor atom can be taken into account by introducing a phase 
shift as 

(f) b (r) = ex PH7^:2/] ~ d %) > ( 37 ) 
where ip is the variational wavefunction in Eq. (J32l) and the corresponding energy is given 
by Eq. ( 1341) . The relevant integrals appearing in the expression of J in Eq. ( |29i) are given 
in Appendix B. 

The Bohr radius ag of the donor electron is about a few nm, whereas the confinement 
length in Eq. (|3~T1) is about 10-20 nm for hu>o = 1 ~ 5meV. Thus the spatial overlap 
between <p a and 0& is expected to be small and correspondingly the exchange integral Vx 
is small, leading to the weak dependence of the singlet-triplet splitting energy J on the 
magnetic field. In this range the Zeeman energy is comparable to or larger than Vx and 
determines the magnetic field at which the singlet-triplet crossing occurs. The regime (b), 
where one electron is bound to the donor and the other electron is delocalized over the QD, 
corresponds to the situation that the harmonic confinement potential is compensated by the 
donor binding energy, namely, 



—mujnd 2 ~ —R 
2 2 



d 



R 

hu 



a B . 



(38) 
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For ZnSe, ob=3.0 nm, R=52.6 meV and thus for hwo=5 (10) meV, the regime (b) cor- 
responds to d ~ 31.6 (15.8) nm. The harmonic confinement length £o = \Jhj (mug) is 
estimated as 9.76 (6.90) nm and dji§ = 3.2 (2.3). The S-T + crossing behaviors are exhib- 
ited in Fig. [5] for ZnSe QDs with g e = 1.2. Figure [5]^a) exhibits the case of huo = 5, 10 
meV and hjj z = 25 meV with a donor at the distance d = 3 £q from the center of the 
QD. A clear S-T + level crossing occurs at an easily available magnetic field strength. With 
increasing fiw z , the wavefunction is squeezed along the z-direction and the spatial overlap 
between wavefunctions and accordingly the exchange Coulomb energy is decreased. In order 
to increase the exchange Coulomb energy, more tighter magnetic confinement of electrons 
is required, pushing the S-T + level crossing to the higher magnetic field for the stronger 
confinement in the z-direction. Figure Mjo) illustrates these behaviors for hwo = 5 meV and 
fkj z =5, 25 and 50 meV with a donor at the distance d = 3 £q. In this case also, the S-T + 
level crossing occurs at around 1 Tesla. 

Finally, we discuss the effect of the spin-orbit (SO) coupling on the S-T-t anticrossing gap 
in the regime (b) of the donor-doped QD. The detailed arguments are developed in Appendix 
A. A unitary transformation is applied in order to eliminate the linear-in-momentum SO 
coupling. Then the most dominant term is the SO induced Zeeman interaction as 
discussed in Sec. II B, which is of the first order in the smallness parameter defined by 
e = £t/^so- Unfortunately, the matrix element (S\H§°\T±) cannot be made to vanish by 
tuning the direction of the magnetic field. Thus the magic angle tuning is not possible 
here. Instead, in order to reduce the matrix element, we have to use the group IV elemental 
semiconductors in which the Dresselhaus SO terms are absent and the magnitude of the 
Rashba SO term might be reduced very much by careful tuning of the strain fields and 
the electric field. If we can achieve it/Xg, £t/\ ~ 10 -4 , the above matrix element would 
be about 10 -8 eV. At the same time, the contribution from the higher order perturbation 
series with respect to H§ is much smaller. The contributions from the renormalized SO 

(3) 

interaction H ren , the cubic-in-momentum SO interaction and their associated higher 
order perturbation terms are estimated in Appendix A and are of the order of 10~ 12 ~ 10~ n 
eV. Thus the S-T anticrossing gap is mainly determined by H§° and can be made to be 
of the order of 10~ 8 eV by using the group IV semiconductors and by tuning the Rashba 
SO coupling constant. On the other hand, the hyperfine coupling energy in the localized 
electron system is about 10~ 7 ~ 10~ 6 eV as discussed in Sec. IV. Thus, in these systems 
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FIG. 5: The energy difference between the triplet (T+) and the singlet (S) states for a donor doped 
ZnSe QD is plotted as a function of the magnetic field. The left panel corresponds to the case (a) 



hujQ = 5, 10 meV and hw z = 25 meV with a donor at the distance d = 3 £o(£q = ^/(mwo)) from 
the center of the QD. The right panel exhibits the case (b) Hujq = 5 meV and huj z = 5, 25, 50 meV 
with a donor at the distance d = 3 £q from the center of the QD. 



our protocol of the nuclear spin quantum memory can be fully effective. 



IV. HYPERFINE INTERACTION AND QUANTUM STATE TRANSFER 

We describe here the protocol for the quantum state transfer (QST) between a pair of 
electrons and a single nucleus. As such examples, there are a single electron charged 28 Si QD 
with a single P atom, a two-electrons charged 28 Si QD with a 29 Si atom, a single-electron 
charged ZnSe QD with a single F atom, a two-electrons charged ZnSe QD with a 77 Se isotope 
atom, and similar structures of the group IV elemental and II- VI compound semiconductors. 

The hyperfine (hf) interaction between a pair of electrons and a nuclear spin is given by 
the contact hf interaction 2 ^ 

Vhf = ^-ge^B 9n^n S« ■ I £(r 4 - R) , (39) 

i=l,2 

where S(I) is the electron (nuclear) spin operator in the dimensionless form, /^s(// n ) the 
Bohr (nuclear) magneton, g e {g n ) the electron (nuclear) g-factor and R denotes the position 
of the nucleus. For the QST protocol we will be interested in the singlet (S) state and the 



1(3 



7 ex P{-7T^} cos -r > ( 4 °) 



triplet (T+) state with the magnetic quantum number m = 1 of an electron pair. Here we 
consider for simplicity the case of two delocalized electrons in a QD containing a single 
isotope atom of the host material with a nuclear spin. The orbital function 4> g (4> e ) is the 
ground (first excited) state with the magnetic quantum number m=0 (1) of the solutions of 
the effective mass equation Eq. (J3]): 

1 [Y . V 2 , 7TZ 

/ 2 T 7TZ 
0e = ^^ r X P {-— -^COS^, (41) 

where i is the harmonic confinement length, d the thickness of QD and the z coordinate lies 
in the range of [— d/2, d/2]. Then the wavefunctions of the singlet S state and the triplet 
T + state are given by 

\S) = ^(^(r^^O/^k) - £(6)a(6)) , (42) 

\T+) = -J=[0 9 ( ri )0 e (r 2 ) - e (n)0 9 (r 2 )] a(&)a(6) , (43) 

where a and /3 are the spin up and down functions and £ denotes the spin coordinate. 

In the vicinity of the S-T + crossing, the hf interaction comes into play, inducing a flip-flop 
process between the electron spin pair and the nuclear spin. The relevant Hamiltonian near 
the crossing point can be derived by calculating the matrix elements of the Vhf. It is to be 
noted that the actual wavefunction of an electron in solids is a product of the Bloch function 
u(r) and the envelope function F(r) in, e.g., Eqs. ( l40l) and (14TT) : 

(r) = F(r) u(r) (44) 
with — / rfr|n(r)| 2 = l, / dr\F(r)\ 2 = 1 , (45) 

v Jv J 

where the Bloch function is normalized in the volume vq of a unit cell and is dimensionless 
and F(r) is normalized in the whole space. Now the relevant Hamiltonian is given as 

V hf = h \T + ) e e (T + \ ® (| |)n nil I " I l)n n(i I) 

+/ii(|T + > e U) ne <S , | n (T|+h.c.) (46) 
with / io = ^(|0 9 (R)| 2 + |0 e (R)| 2 ) , (47) 

ht = ^:(R)0,(R) , (48) 

A = ^-g e [i B gnfJ>n\u(K)\ 2 , (49) 
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where R denotes the site of the nuclear spin and the suffix e(n) attached to the ket and bra 
vectors indicates the electron (nucleus). 

For the single electron charged QD doped with a donor atom, we can distinguish three 
regimes according to the donor position as discussed in Sec. III. In the regime (c) where 
two electrons are delocalized throughout the QD, the hf coupling Hamiltonian is exactly the 
same as obtained above. In the regime (a) where two electrons are strongly bound by the 
donor atom, the same formulation as above can be used only by changing the meaning of 
the basis functions <fi g and <fi e . Namely, <p g and <p e are the donor-bound ground and excited 
states, respectively. Then hi vanishes because <j) e is usually an odd-parity excited state which 
has a vanishing amplitude at the origin, i.e., at the donor nucleus. On the other hand, in 
the regime (b) one electron is delocalized within the QD and the other electron is strongly 
bound by the donor atom. In this case the singlet and triplet states are given by 

\S) = ^(rO^ra) + &(ri)0«(r a ))(a(£i)/3(6) - , (50) 

\T+) = -L[0 a ( ri )0 6 (r 2 ) - 6 ( ri )0a(r 2 )] <i)a(6) , (51) 

where <f) a is the ground state in the QD confinement potential and <f>b is the ground state of 
the donor-bound electron. Then we find 

^ = ^(IMR)| 2 -|0a(R)| 2 ). (52) 

Noting that in the regime (b) the donor atom is located off the center of the QD potential by 
a few times the confinement length £, we can neglect the second term in Eq. (1521) . The first 
term can be large because ^ is a strongly localized function with a typical spatial extent of 
a few nm. For example, in a P-doped Si QD h\ would be 41.4 MHz corresponding to the 
time constant of 24 ns, because v4|0b(R)| 2 is known to be 117 MHz from the experiments 
on Si:P samples 2 ^. Consequently, the regime (b) of the single-electron charged QD doped 
with a donor atom is most favorable to realize a strong hf coupling and to achieve the 
electron-nuclear spin QST. 

Now we shall discuss the feasibility of the electron-nuclear spin QST for the two-electrons 
system. For the two-electrons charged ZnSe QD doped with an isotope 77 Se atom having 
the nuclear spin(/=l/2), the hf coupling constant would be about A = 3 ~ 4 x 10 -8 eV 
(nm 3 ), if \u(R)\ 2 = 100 is assumed for the Bloch function referring to the value 186 in the 
case of Si^£. For a typical QD with £ ~ 20 nm and d ~ 10 nm, the magnitude of the hf 
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interaction hi is estimated as 



h x = 2 ~ 3 x (|R|/£)exp[-|R| 2 /£ 2 ] x 10~ 12 eV ^ 2 ~ 3 x 1(T 12 eV , (53) 

where the nuclear spin is assumed to be located at the midpoint z = in the z-direction 
and the factor related to the lateral part of wavefunctions is simply assumed to be 1. The 
corresponding electron-nuclear spin QST time is 1~ 2 ms. The singlet (S)-triplet (T) anti- 
crossing gap due to the SO interaction should be much smaller than the hf coupling energy. 
As shown in Appendix A, the S-T anticrossing gap is determined by the matrix element 
(S\VgQ |T±), where V^J is the effective SO interaction after a unitary transformation to 
eliminate the original SO interaction. The most dominant term is the SO induced Zeeman 
interaction H§° and this term can be made to vanish in the first order by the magic angle 
tuning of the direction of the magnetic field. But the actual value of the S-T anticrossing 
gap is determined by the higher order perturbation terms concerning H§° and contributions 
from the renormalized SO interaction H ren . The typical magnitude of the S-T anticrossing 
gap is about 1CT 12 ~ 1CT 11 eV as shown in Appendix A. Thus the S-T anticrossing gap 
is of the same order of magnitude as the hf interaction energy. On the other hand, in Si, 
the Dresselhaus SO terms are absent due to the inversion symmetry of the crystal lattice. 
According to the experimental report on a SiGe/Si/SiGe quantum well^I, the Rashba SO 
coupling constant Han is about 0.55 xlO -4 eV A. In this case also, the S-T anticrossing 
gap is estimated to be about 10~ 12 ~ 10~ n eV. The hf coupling constant for the 28 Si QD 
with a 29 Si(/=l/2) isotope atom is calculated as A = 5.06 x 10^ 8 eV (nm 3 ), employing 
|w(.R)| 2 = 186^ and hi is estimated as 

h x = 4.0 x 10~ 12 eV (54) 

for I = 20 nm and d — 10 nm. Thus, in Si and ZnSe, the S-T anticrossing gap is of the 
same order of magnitude as the hf coupling energy. Consequently, our protocol for the 
electron-nuclear spin QST will not be effective. 

On the other hand, in the single-electron charged QD with a donor impurity having 
the nuclear spin, the stronger hf coupling and the faster electron-nuclear spin QST can be 
expected because of the highly localized nature of the donor wavefunction. In the single- 
electron charged Si QD doped with a P donor, two electrons composed of the donor electron 
and an externally added electron play the role of a qubit. As discussed above, the hf 
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coupling is the strongest for the regime (b) of Fig. HI in which one electron is delocalized 
within the QD and the other electron is strongly localized around the donor atom, and 
the magnitude would be 41.4 MHz. Concerning the spin-orbit interaction, the Dresselhaus 
terms are absent and the Rashba term gives rise to the S-T anticrossing gap of the order of 
10 -8 eV, as shown in Appendix A. For the single-electron charged ZnSe QD doped with a F 
donor, the situation is the same, although the S-T anticrossing gap would be larger. For this 
system the hf coupling constant would be a few tens of MHz according to the estimation 
in Appendix D. Thus, the single-electron charged QDs of the group IV elemental and II- VI 
semiconductors doped with a donor impurity is very favorable to realize the electron-nuclear 
spin QST and the nuclear spin quantum memory. 

In order to realize QST, a flip-flop type of interaction is essential and in a QD occupied 
by a pair of electrons this can be realized at the S-T + crossing. When the nuclear spin is 
initialized in the j state, the system evolves as 

1 (,(t = 0) = (a\S) e + b\T + ) e )®\l) n , (55) 
ip(t) =a\S) e \ l) n + b(cos-^ \T + ) e \ Dn-isin^ \S) e \ f) n ) , (56) 

V'(< = |^) = |5>«®(oU>n-i6|T>»), (5T) 

where a and b are arbitrary constants normalized as |a| 2 + |6| 2 = 1. Thus the quantum state 
of the electron pair is transferred to the nuclear spin. After that by tuning the magnetic 
field strength off the S-T + crossing point, the hf interaction is effectively switched off and 
the nuclear spin memory can be preserved. In the retrieval process, we prepare the singlet 
state of an electron pair in the QD and tune the magnetic field strength just at the S-T + 
crossing point. Then by waiting for a time nh/(2hi), the state evolves as 

^(0 = \S) e ® (a| I) n - ib\ T>„) - ipit' + = (a\S) e - b\T + ) e ) ®\[) n . (58) 

This state is not exactly the original state but the sign change of b can be remedied by the 
optical STIRAP process which rotates the pseudospin spanned by the \S) and |T + ) states^. 
If we want to recover the original state only by the hf interaction, we have to wait for a 
time 3irh/(2hi). In the above we have neglected the nuclear Zeeman energy. When this 
Zeeman energy is taken into account, another phase factor is attached to b in Eq. (1561) . This 
phase factor, however, can be cancelled by inverting the direction of the magnetic field in 
the retrieval stage. 
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In this scheme, a knowledge of the hf interaction constant hi is necessary to achieve QST. 
Now we discuss how the hf interaction constant hi can be determined. This is a sample 
specific information and can be fixed by measuring the probability of the spin state mixing 
between the singlet and triplet states of two electrons as a function of the hf interaction 
period, in a similar way to the experiment carried out by Petta et al.— . One has to initialize 
the two electrons in the singlet state, then to adjust the magnetic field strength to the 
singlet-triplet crossing point and to wait for a period r. Next, tuning away from the singlet- 
triplet crossing point, one measures the probability to find the electrons in the singlet state. 
The measurement is repeated at a long interval during which the nuclear spin is randomized. 
The probability to obtain the outcome of the singlet state should follow 



P(t) 



1 



2( hiT, 

cos 



(59) 



2L v h 

from which the magnitude of the hf interaction constant hi can be estimated. 

Nuclear spin initialization can as well be realized through the electron spin measurements. 
Assume that a pair of electrons is initialized in the spin singlet state, whereas the nuclear 
spin is in an arbitrary mixed state, then the system is tuned to the S-T + crossing point for 
a time r. The time evolution of the system is given by 

P (t = 0) = \S) e e (S\ <g> f PT I |)n „<T \ +Pi I On nU l) > (60) 



p(r) = P] Mr)) Mr) I + Pl \S) e e (S\ <g> | |) n „(| | (61) 

with V(r) =cos^ \S) e \ t) n -isin^ \T + ) e \ |) n . (62) 
n n 

At the next step an electron spin measurement is carried out in the singlet-triplet basis. The 
probability Pg to have the outcome of the singlet state and the density matrix p s after the 
measurement are given by 

h t 

P.s=Pi+ P T cos2 -jr- > ( 63 ) 



PS = ^\S)e e(S\ ® ( H I |) n n (| I +PT COS 2 ^ | |}n n(T |) • (64) 



In the same way, for the outcome of the triplet state the corresponding quantities are given 
by 

P T = p T sin 2 — , p T = \T+) e e (T+\ ® | |) n n {i | . (65) 

If one knows the hf interaction constant, then setting r = irh/ (2hi) guarantees that the 
nuclear spin is initialized in the j state. Otherwise, i.e., when one doesn't know the exact 
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FIG. 6: Schematic energy level diagram around the S — T + crossing point. Dashed lines (|T+) e | |) n 
and \S) e \ |) n ) represent the original energy levels without the hf coupling, whereas solid lines 
(\E + ) and \E_)) depict the energy levels including the hf coupling. A denotes the energy difference 
between |7+) e | [) n and \S) e \ j) n . The adiabatic state transfer is achieved by sweeping the magnetic 
field from B{ to Bf and backward. 

hf interaction constant, the initialization goes as follows: if the triplet state is detected, the 
nuclear spin state is initialized in the J, state deterministically. If the singlet state is detected, 
in the conditional state ps(Eq. (1641) ). the weight of the f nuclear spin state is decreased 
with respect to the initial state p(t = 0)(Eq. (I60p). Thus by continuing the electron spin 
measurement, one can eventually purify the nuclear spin state into the j spin state. 

V. ADIABATIC QUANTUM STATE TRANSFER 

In the last Section, the scheme of the quantum state transfer between two electrons and 
a nuclear spin is presented under the assumption that the hf coupling constant is known. 
However, there is an alternative scheme to do the same quantum state transfer without the 
knowledge of the hf coupling constant. This is based on the adiabatic state control by the 
magnetic field tuning through the S-T + crossing point, which is schematically illustrated in 
Fig. m 
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The relevant Hamiltonian is given by 

H=Z Y ^ ^ (I t>» n{\ I + I l>n n(i |) + y|T + ) e e (T + \\ |) n n (| | 

(\T + )e\ i)n e(S\ „(? | + \S) e \ T>„ e (T + 1 n (| |) , (66) 

where A is the energy splitting between the singlet and triplet states, hi the hf interaction 
constant derived in Eq. (1481) and the nuclear Zeeman energy is omitted. The electron spin 
states \T + ) e and \S) e are coupled by the hyperfine interaction, resulting in the eigenstates: 

\E ± ) = JL[-/i 1 |r + ) e | l} n + (A/2 T E )\S) e \ |)J (67) 



with A± = hi + (A/2 =f £ ) 2 (68) 



and the energies 



E ± = ±E , E = JAy4 + hj. (69) 



Initially an arbitrary electron state is prepared with the nuclear spin in the | |) n state. 
Then the magnetic field is sweeped to the right from B — B t at t — ti to B — Bf at t — tf, 
as shown in Fig. 5 adiabatically and the parameter A is tuned from an initially positive 
value Aj to a negative value A/, where |Ajj| ^> h\ is assumed. At B — B t or A = A; the 
state is composed mostly of \T + ) e \ |) n , whereas at B = Bf or A = Af the state \E + ) 
is composed mostly of \S) e \ t}«- In the course of this sweeping, the \T + ) e \ |) n state starts 
from the \E + ) branch in Fig. 5 and continues to be on the same branch but an extra phase 
factor is acquired. The time evolution of the wavefunction is given by 

l^o) = [a\S) e + b\T + ) e }\ l) n - = \S) e [a\ i) n e^+b\ t)„e^ 2 ] (70) 
1 f tf 1 f Bf 1 



U JBi | dt 



with ct> l = - dtA(B) = -l dB-^A(B), (71) 



02 = ~\ f " dt yjA*(B)/4 + hl = -\ I*' dB-^ y/A*(B)/4 + h* , (72) 

where A is a function of the magnetic field B, B is a function of the time t and hi is assumed 
to be constant. As a consequence, an additional phase difference 02 — <Pi is introduced in the 
superposition state of the nuclear spin. However, this additional phase can be eliminated 
at the retrieval stage. A radio frequency 7r pulse is applied on the state \ipi), swapping the 
nuclear spin states | |) n «-> | |) n and the state becomes 

|^ 2 > = exp^vrM) = \S) e [ a\ ]) n e^ + b\ [) n e^\ , (73) 
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where I x is the x-component of the nuclear spin operator (|/| = 1/2). 

Then the adiabatic sweeping is reversed from Bf at t — to Bi at t — t'^ with the same 
time profile. In this case, the state \S) e \ |)„ transfers adiabatically to the state \T + ) e \ |) n , 
although another extra phase factor is acquired: 

02 = ~\ dt y/A*(B)/4 + hl = -\ dB-^ y/A*(B)/4 + h* = . (74) 

Jt'i J Bi \ dt\ 

In the same way, the state \S) e \ |) n acquires the phase shift 

= in r dtA(B) = h r a(b) = 01 • (?5) 



I d.B | 
dt 



Consequently, the wavefunction becomes 

|^ 2 > - a\T + ) e \ I) n e^+^ + b\S) e \ I) n e^ + ^ = [a\T + ) e + b\S) e ]\ [) n e^ + ^ . (76) 

This state is not just the original state |^o) in Eq. (1701) but with a swap between the spin 
singlet and triplet states. However, the above protocol can be repeated once more to obtain 
the original electron spin state. Alternatively, we can employ an optical STIRAP (stimulated 
Raman adiabatic passage) method to perform the rotation within the pseudospin subspace 
composed of the states |S') e and \T + )J&. Then the original state is recovered. 

Finally, the condition for the adiabatic state transfer will be discussed. When the mag- 
netic field is sweeped, the parameter A in the Hamiltonian ( 1661) is time-dependent. The 
wavefunction in the subspace spanned by \S) e \ f)n and \T + ) e \ |) n can be expanded as 

\m) = a(t)\E + ) + b(t)\E.) , (77) 

where the eigenstates \E + ) and \E_) are also time-dependent through the time-dependence 
of the parameter A. By noting that 

J t \E + ) = ^j&M - A / 2 )l^> = C +W \ E ~) • ( 78 ) 

J t \E-) = + A / 2 )l^) = c " W 1^) > ( 79 ) 

where c + (t) and C-(t) are introduced for simplicity, we have the Schrodinger equations: 

j t a(t) = ~E + (t) a(t) + c_(0 b(t) , (80) 
jb{t) = - l -E.{t)b{t) + c + {t)a{t). (81) 
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When the second terms on the right hand side of the above equations can be neglected, 
the adiabatic time-evolution is guaranteed. This condition will be examined near the S-T + 
crossing point where A ~ and the adiabaticity condition is most stringent because the \E±) 
branches are close to each other. Around the S-T + crossing point, E + = hi , £L = —hi 
and c + (t) and c_(£) are estimated as 

c + (t) = -c„(t)- A = CQ . (82) 

Then the time evolution under the initial condition a(0) = 1 and b(t) = is given by 

a(t) = cos Qt — i-^- sin Qt , b(t) = sin Qt (83) 



with ui = -± , ft = \JujI + cI . (84) 

The amplitude b(t) should be small enough to guarantee the adiabaticity, namely, 

. A , , hi , , : , , dB dA . Ah 2 , 

77- < ^ — ► A = — — < — i. 85 
Ahi h dt dB h 

This condition restricts the speed of the magnetic field sweeping for the adiabatic state 
transfer. 



VI. QUANTUM STATE TRANSFER BETWEEN A SINGLE ELECTRON AND 
A SINGLE NUCLEAR SPIN 

So far we have considered the qubit composed of a pair of electrons. However, the 
single electron qubit is more fundamental as a building block for the quantum information 
processing. Here we consider a system composed of a single electron as a qubit and a nuclear 
spin as a quantum memory whose examples will be discussed later. It is possible to devise 
a scheme for the quantum state transfer between a single electron spin and a nuclear spin. 
At zero magnetic field the hyperfine interaction between a single electron spin and a nuclear 
spin is given by 

H = A I • S , A = ^g efiB9nfin |F(R)w(R)| 2 , (86) 

where I(S) denotes the spin-1/2 nuclear (electronic) spin, u(r)(F(r)) is the Bloch (envelope) 
function of the electron and R is the position vector of the nucleus. In the hf interaction 
the dipolar term vanishes, since the dynamics takes place within the electronic ground state 
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orbital and the Fermi contact interaction is dominant. A magnetic field is applied in order 
to switch off the hf coupling so that the large difference in the Zeeman energy splitting 
between an electron and a nucleus prohibits the flip-flop transitions between them. As 
shown in the Appendix C, under the Hamiltonian in Eq. (|86j) . the time-evolution of the 
electron-nucleus coupled system proceeds as follows: 

p (t = 0)= \tp) e e(lp\ ®Pn ►/£>(* = nh/A) = Pe ® \<lp) n n (^\ (87) 

with |^> = a|T>+m>, (88) 

where a and f3 are arbitrary complex constants normalized as | or | 2 + \(3\ 2 = 1, the suffix 
e(n) indicates the electron (nucleus), p n is a density matrix representing an arbitrary mixed 
state of the nucleus and p e is the same mixed state for the electron. Thus the quantum state 
transfer is possible for an arbitrary mixed state of the nuclear spin and the initialization of 
the nuclear spin is not necessary. 

In the first stage of the QST, after the hf interaction of duration r = nh/A the spin 
states are swapped between the electron and the nucleus. Then a magnetic field is applied 
to switch off the hf coupling and to preserve the nuclear spin memory. After some interval 
within the nuclear spin coherence time, the magnetic field is turned off and the hf interaction 
is again switched on for another time period of r = nth/ A. Then the initial quantum state is 
retrieved back in the electron spin. Actually, the magnetic field applied during the storage 
period introduces an extra phase factor due to the nuclear Zeeman energy splitting. But 
this phase can be cancelled by applying a magnetic field after retrieving the quantum state 
in the electron spin. 

This scheme of QST is relevant for a single electron trapped in a single QD and also for 
a neutral donor in the bulk semiconductor which is composed of isotope atoms without the 
nuclear spin. For instance, the Si:P system has been extensively studied and has already 
been proposed as a qubit for the scalable quantum computer— 31 P has a nuclear spin 
/ = 1/2 with the hf interaction constant of A = 117 MHz. Thus the QST can be realized 
on a time scale about 10 ns, which is much shorter than the electron coherence time which 
can extend to about 60 ms^^. Another candidate might be a II- VI semiconductor because 
in this material the abundance of isotope atoms with the nuclear spin is only a few percent 
and the isotope purification may be possible. For example, ZnSe has the natural abundance 
of 4.1(7.6)% for the 67 Zn ( 77 Se) isotope atom with the nuclear spin I = 5/2(1/2). Thus the 
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isotopically purified and nuclear spin-free ZnSe crystal doped with a fluorine atom, namely 
ZnSe:F, would be a good example^. So far, we are not aware of any NMR data for the 
ZnSe:F system. But we can estimate the hyperfine coupling constant for the ZnSe:F by 
comparison with the case of Si:P, as shown in Appendix D. According to that, we infer that 
A(ZnSe:P) ~ 78.6 MHz, which indicates a smaller hf coupling compared with the Si:P case. 

The important feature of these localized electron system is the optical interface. Because 
of the localized nature the optical transition is possible even in the indirect-gap materials 
like Si. Actually the photoluminescence from Si:P centers has been extensively studied^ 1 ^. 
The photoluminescence occurs through the neutral donor-bound exciton (D°X) state which 
is composed of two electrons and one hole. This photoluminescence can be utilized to locate 
and address each localized center. In view of the recent progress in manipulating optically 
a single electron spin trapped in a QE>2L2£, we can expect that the donor electron spin can 
also be controlled optically via the A-type transition through the donor-bound exciton state. 
Another important feature of the localized electron system is the uniformity of the system. 
For example, any Si:P center has the same optical transitions like atoms because the nearby 
atomic configuration is the same for any center. Thus a laser light of the same wavelength 
can be used for the optical initialization, manipulation and measurement of the electron spin 
in any localized center. This feature is advantageous in constructing a scalable quantum 
network from an array of these localized centers through the optical interconnection. 



VII. SUMMARY AND DISCUSSION 



We have proposed a few new protocols for the nuclear spin quantum memory in the 
isotopically purified group IV elemental and II- VI compound semiconductors, where the 
number of atoms with the nuclear spin is reduced to only one and thus the decoherence due 
to the nuclear dipole-dipole interaction can be avoided. We have studied various cases where 
the qubit is defined by the singlet and one of the triplet states of a two-electron system or 
simply by a single electron spin. For the two-electron system, we studied the case where two 
electrons are loaded on a donor-free QD with an isotope atom having the nuclear spin and 
also the case where one electron is additionally loaded on a QD doped with a single donor 
atom. In the latter case the behavior of the singlet-triplet crossing depends sensitively on 
the donor position, reflecting the hybridization between the donor-bound localized orbitals 
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TABLE I: Properties related to the quantum state transfer (QST) for the two-electron and single- 
electron qubits. Schematic configurations of electrons and an impurity for cases (a) to (f) are 
illustrated in Fig. 7. 





two-electron qubit 


single-electron qubit 


(a) Si QD 


(b) ZnSe QD 


(c) Si:P QD 


(d) Si:P 


(e) ZnSe:F 


(f) ZnSe QD 


Nuclear spin 
initialization 


necessary 


necessary 


necessary 


unnecessary 


unnecessary 


unnecessary 


Electron- nuclear 
QST time 


~ 1 ms 


~ 1 ms 


~ 20 ns 


~10 ns 


~ 10 ns 


~ 1 ms 


Spin-orbit vs. hy- 
perfine coupling 


Vso ~ Vhf 


Vso ~ Vhf 


Vso <V hf 


irrelevant 


irrelevant 


irrelevant 


S-T-t- crossing 


unfavorable 


unfavorable 


favorable 


irrelevant 


irrelevant 


irrelevant 


Optical Interface 


hard 


possible 


good 


good 


good 


possible 



(a) 28 Si QD (b) ZnSe QD (c) 28 Si QD 




(d) 2S S j bu | k ( e ) ZnSe bulk (0 ZnSe QD 




FIG. 7: Schematic configurations of electrons and an impurity for the electron-nuclear spin QST 
which are compared in Table I. The solid circle (square) represents an electron (impurity). The 
cases of (a), (b) and (c) correspond to the two-electron qubit, whereas the cases of (d), (e) and (f) 
to the single-electron qubit. 
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and the delocalized QD orbitals. Our protocol requires that the hyperfine coupling energy 
should be much larger than the singlet-triplet anticrossing gap. This requirement is satisfied 
favorably in the case of a donor doped QD because the large hyperfine coupling energy 
can be achieved only for strongly localized electrons. In the two-electron system the 100% 
nuclear spin initialization is possible based on the spin state measurement of electrons and 
can be utilized in the QST between the electronic and nuclear spins. Our protocol of the 
nuclear spin quantum memory can be achieved by the electrical control even without tuning 
the magnetic field and complements the advent of the electrical gate tuning of the singlet 
and triplet states of a pair of electrons^. 

For the single electron system, we considered the case of a single donor electron coupled 
with the nuclear spin of the donor atom in the bulk crystal or in a QD. Here the QD is not 
essential but is convenient only for addressing each electron locally, because the donor elec- 
tron is much more tightly bound to the donor atom compared with the delocalized electron 
in a QD. The most important feature in the single electron system is that the nuclear spin 
initialization is not required in the QST between the electronic and nuclear spins. Typi- 
cal examples of the single electron system are Si:P and ZnSe:F. These systems have the hf 
coupling constant about a hundred MHz and the QST time between the electron and the 
nucleus is about 10 ns. This feature is favorable in view of the long coherence times of the 
electron spin and the nuclear spin in the group IV elemental and II- VI compound semicon- 
ductors. Another example is the case of a single delocalized electron in a QD containing a 
single isotope atom with the nuclear spin of the host material, e.g., an isotopically purified 
28 Si QD with a 29 Si atom and an isotopically purified ZnSe QD with a 77 Se atom. Although 
the hf coupling constant for some cases is not well known, the constant can be inferred by 
assuming that \u(R)\ 2 ~ 100 in Eq. ()49p . According to this rough estimate, the hf coupling 
constant for a delocalized electron in a QD with a typical size of 20 nm is about a few 
times 10~ 12 eV and the QST time is about 1 ms. Thus the delocalized electron in a QD 
is not favorable compared with the donor-bound localized electron for the application to 
the nuclear spin quantum memory system. Furthermore, the donor-bound localized electron 
has, in general, a good optical interface and the optical initialization, manipulation and 
measurement of the electron spin would be possible through the donor-bound exciton state, 
although the experimental demonstration is yet to be challenged. 

In Table HJ merits and demerits of the two-electron qubit and single-electron qubit are 
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compared with respect to the time required for the QST between the electron and nuclear 
spins, the necessity of the nuclear spin initialization in the QST, the relative magnitude be- 
tween the spin-orbit coupling and the hf interaction which determines the feasibility of our 
QST protocol, the possibility of the singlet-triplet state crossing for the two-electron qubit, 
and the feasibility of the optical interface. From this Table, we see that the donor-bound 
localized electron system is favorable with respect to the fast QST between the electron and 
nuclear spins, the irrelevance of the nuclear spin initialization in the QST and the feasibility 
of ultrafast optical manipulation and measurement of the electron spin. Furthermore, the 
localized electron system is homogeneous like atoms in the sense that the energy level struc- 
tures and the associated optical transitions are the same for any localized center because 
of the characteristic nearby atomic configuration of a specific localized center. Based on 
these features, we envisage that an array of the donor-bound localized electrons provide an 
excellent set of qubits and the gate operation or entanglement transfer between any two 
localized qubits can be carried out through optical channels, as demonstrated recently using 
atomic qubits^, and consequently the quantum network can be established. 
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APPENDIX A: SPIN-ORBIT INTERACTION FOR A PAIR OF ELECTRONS 



The spin-orbit (SO) interaction for the conduction band electron in the linear approxi- 
mation with respect to the momentum operator is given by 

V S o = a R (a x p y - a y p x ) + a D (a x p x - a y p y ) (Al) 

in the two-dimensional limit, where the first term is the Rashba term due to the structural 
inversion asymmetry and the second term is the Dresselhaus term arising from the bulk 
inversion asymmetry. In the case of a pair of electrons it is convenient to rewrite the Hamil- 
tonian in terms of the center-of-mass coordinate R = (X, Y, Z) and the relative coordinate 
r = (x, y, z), defined by 

R=-(r 1 + r a ) , r = ri -r 2 , (A2) 

where the subscript 1(2) refers to the first (second) electron. Then we have 

d d d d 1 d d 

and it is natural to introduce the momentum operators corresponding to the center-of-mass 
coordinate and the relative coordinate by 

n = P i + p 2 , P = i( P i-p 2 ) . (A4) 

Under a uniform magnetic field the momentum operators are modified into the gauge- 
invariant form: 

Pi = -ifiS- + -A(ri) = ~ih^- + ^Bx ri (i = 1, 2) , (A5) 

OTi C OTi 2c 

n = -ih— |-BxR, (A6) 

Oil C 

P = + x r , (A7) 

or 4c 

where the symmetric gauge (A(r) = (B x r)/2) is employed. Accordingly, the orbital part 
of the Hamiltonian is rewritten as 

h = -L(n2 + n; + n*) + u cm (x, y, z) + -( P 2 X + v \ + v \) + u rd (x, y,z) + ^- (As) 

where U cm (U re i) are the circularly symmetric confinement potential for the center-of-mass 
(relative) coordinates and it is to be noted that the mass of the center-of-mass (relative) 
coordinate is 2m(m/2). 
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Now the spin-orbit interaction in Eq. (IAip and the Zeeman energy Hz for two electrons 
can be written as 

Vso = ^an^Jly - Y, y Il x ) + a R (a x p y - a y p x ) 

+ ^a D (E x Il x - Zylly) + a D (a x p x - a y p y ) , (A9) 

H z =-^sB-S (A10) 
with Si = an + a-n , <7, = a u - a 2i (i = x, y, z) , (All) 

where II and p are defined in Eq. (1A4I) . Now we examine the effect of this spin-orbit inter- 
action on the singlet-triplet level crossing. Since the Hamiltonian in Eq. (IA8I) is separated 
into the center-of-mass and the relative coordinates, the singlet and triplet eigenstates can 
be written as 

\S) = ^ cm (R)^(r)i=( a (a)/3(6) - /?(6M6)) , (A12) 
\T + (T_)) = g cm (K)e rel (v) O3(£i)0(&)) , (A13) 

|T ) = &cm (R)e^(r)-^=(a(6)/?(6) + /?(£iM&)) , (A14) 

where g C m{.9rei) is the ground state of the center-of-mass (relative) motion, e re / is an odd- 
parity excited state of the relative motion, a(£) and /3(£) are the spin up and down functions 
and £ denotes the spin coordinate. The total spin operators S x and E y do not change the 
magnitude of the total spin but they and the difference spin operators a x and a y change the 
magnetic quantum number only by ±1. Thus the following matrix elements vanish: 

(S\V so \S) = (T ± \V so \T ± ) = (T ± \V so \T T ) = (S\V SO \T ) = (T \V SO \T ) = . (A15) 

In the matrix element (S\ Vso\T±), S x and E y do not contribute because they do not change 
the magnitude of the total spin and thus only the a x and a y contribute. After the spin part 
is calculated, we find 

(S\V so \T±) = (g C mgrei\i\/2a R (p x ± ip y ) =f V2a D (p x =f ip y )\g cm e r ei)- (A16) 

Then, using the relations: 

[x ' H] = ^T2 v " [v ' H] = ^' (A17) 
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we have 

Tfi 

(g cm grel\Px ± iPy\gcm.e r el) = ^ {9cm9rel \ ± H]\g cm e re i) 

= 2 ^ (prefix d= gy|e re ;) = — (g re /|re ^|e rei } (A18) 
with A ST = E Z {S) - E Z {T ± ) , (A19) 

where Ez{\) denotes the Zeeman energy of the A state. Finally, the matrix element is 
calculated as 

(S\V so \T±) = — ^M^lre^le^) ±ia D {g rel \re^\e rel )} . (A20) 

This matrix element is finite in general because the S-T-t crossing occurs at a finite magnetic 
field. 

However, we have to consider the higher order perturbation terms with respect to the spin- 
orbit interaction Vso- For that purpose it is convenient to apply a unitary transformation 
to the single electron Hamiltonian: 

H = ^—p 2 + a R (a x p y - a y p x ) + a D (a x p x - a y p y ) + ]-g e fi B B ■ cr , (A21) 

where p is the kinetic momentum vector defined by Eq. (1A5I) . By rotating the coordinate 
system by an angle 7r/4 in the xy plane, namely, by introducing a new coordinate system 
(£,?7, z) defined by = (1, 1, 0)/\/2, = (— l,l,0)/v2, e 2 = (0,0,1), the relevant vector 
components are transformed as 





(A22) 

and the same relations hold for the magnetic field B and the vector potential A. Then the 
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above Hamiltonian is rewritten as 



H = — p 2 - Pz<r v (a D + a R ) + p v a^(a R - a D ) 

+ ^fi , e /M 5 £ cr 5 + B v cr v + B z a z ) 
1 / h 



2m 



— \Pv + -^) + ^p 



(A23) 
(A24) 

m(aR + an) ' n rn(aR — Od) ' ^ 
where the last constant term in Eq. (1A24|) will be omitted hereafter. Now, in order to 
eliminate the spin-orbit coupling, we introduce the unitary transformation^: 



-^defJ'BiB^ + B v g v + B z a z ) - m(a 2 D + a|) 



with A 



h 



A,, 



h 



H = H U with C/ = exp[i— av, — i— cr^] . 



A r 



Using the formula: 



• 2 -3 

exp[iS] if exp[-iS] = if + f [5, if] + ^y[5, [5, if]] + |j [5, [S, [S, H}}} 



(A26) 



(A27) 



we calculate the terms up to the second order in S and obtain 
if 



1 1 

T^P 2 + ^9e^B(B ( a^ + B v a v + B z a z 



+9e^B 

h 



V 



A, 



-a c + ^a v )B z -[ —B 6 + —B 



V 



A, 



Ar? V 



-L z a z + ^(-2^ + vit,Pt}) + ^(^Pt ~ t{0,Pr,}) 

A 5 A V 



(A28) 



where {A, f?} = AB + f?A, the first line represents the single electron Hamiltonian under a 
magnetic field without the SO coupling, the second line the SO induced Zeeman interaction 
and the third line stands for the renormalized SO interaction. The energy level diagram 
of H is totally the same as that of the original Hamiltonian H . Thus we can discuss 
the S-T-t crossing behavior based on the transformed Hamiltonian H . It is important to 
note that owing to the unitary transformation the second and third lines of Eq. (1A28|) 
contain a smallness parameter defined by e = it/^M which is typically about 1CT 3 , where 
t t denotes the lateral extent of the electron wavefunction. Thus we can treat these terms 
perturbationally. 



First of all, we study the second line of Eq. ( 1A28I) . which will be denoted by H z 



so. 



uSO _ _ ,. 

U Z — 9e^B 



B, 



B f 



t-B v 

A„ 



(J, 



(A29) 
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For two electrons this part can be rewritten as 

9e^B 



1 



A 



A, 



3 



-t—Btj I 

A^ 



A, 



-CT S 



(A30) 



where R{(£) and Rti(t)) denote the center-of-mass (relative) coordinates in Eq. (1A2I) and 
S^(ct^) and H^a^) are defined by Eq. ( 1A11I) . Then the matrix element between the S and 
T± states is given by 

(T ± \H S Z °\S) = \g e MT±\ + B z - (J-Bz + ^B^ a z \S) . (A31) 

The direction of the magnetic field will be taken as (sin 9 cos (p, sin 9 sin ip, cos 9) in the original 
coordinate system (x, y, z) and is assumed as the direction of the spin quantization. Then 
the spin part of the matrix element is calculated as 



{T± | <?£ | S) = ^V^cos^cosy?- — zv^sin^- , 
(Tilcr^S 1 ) = ^v^cos^sin^. + z\/2cos<£_ , 
(T±\a g \S) = iv^sinfl , 



(A32) 
(A33) 
(A34) 



where y_ = tp — tt/4 is the azimuth of the magnetic field in the (£, rj, z) coordinate system. 
The excited orbital state associated with \T±) can be approximated by 



e r ei(t„v) oc (^ ± irf) g re i{C,v) 



(A35) 



where g re i is the ground state orbital of the relative coordinates. Hereafter, we consider the 
case of £ — ir) in the above equation. Then we have 



{e re j|£|flw) = no , {e re i\v\9rei} = i r 10 , 
where r 10 is a real constant. Consequently, the matrix element is calculated as^ 



(T±\H* u \S)=^-j=g f i B Br ia 



1 cos 9 s 



A, 



± 



cos 9 



sin (p_ 



(A36) 



(A37) 



This result suggests that the matrix element can be made to vanish by tuning the direction 
of the magnetic field. Depending on the relative magnitude between A^ and A^, one of the 
factors 

1 , COS0 . 1 , COsfl /»nn\ 

± — — and — ± — — (A38) 



A 



A, 



A, 
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can be made to vanish by choosing 6 appropriately. Then the factor sin or cos </?_ asso- 
ciated with the non-zero prefactor can be made to be zero by choosing appropriately. 
Usually one of the SO coupling constants ao and cir is much larger than the other and 
thus lA^/AJ = 1. This means that the appropriate angle 9 is nearly or n, indicating the 
^-directed magnetic field. 

The higher order contribution is present. For example, the second order term is given by 

^ {T±\H§ \m}(m\H§ \S} fk . 

£ E { S)-eU (A39) 

where \m) denotes appropriate intermediate states. The typical magnitude of H§° is esti- 
mated for g e ~ 1 and B ~ IT as 

H§° ~ ^ e/ i BJ B A. „ 5 ■ 10" 8 eV , (A40) 

where A^o = A^ = A^ ~ 10/im and £j ~ lOnm are assumed and then the magnitude of the 
above second order term is about a few times 10~ 12 eV because E(S) — E(m) ~ lmeV. 
Now we discuss the renormalized SO interaction given by the third line of Eq. ( 1A28I) : 

The first term L z a z in the parenthesis for two electrons can be written as 



(A41) 



L\ z (Ji z + L 2z a 2z = 2^ Llz + ^2z)£ 2 + 2^ 12 ~ -^Vz , (A42) 

where the subscript 1(2) denotes the first (second) electron. Since the total spin operator 
S z conserves the magnitude of the total spin, we have 

(T±\L lz a lz + L 2z a 2z \S) = -{T ± \{L lz - L 2z )a z \S) 



a 



\S) . (A43) 



This matrix element vanishes when the magnetic field is applied in the z direction, because 
the a z operator does not change the magnetic quantum number. Even in the case of a tilted 
magnetic field, this matrix element vanishes when the orbital functions are given by 

^t±(% R rn £> V) = e cm (% Rn)9rei(£, v) or 9cm{R^ R v )e rd {i, Tj) , (A44) 
i[>s(Rt, Rr» f , rj) = 9cm{Rt, Rr,)g re i(£, v) , (A45) 
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where g is the even parity ground state orbital and e is the odd parity excited state orbital. 
Thus in the first order perturbation the L z a z term can be neglected. Now the higher order 
perturbation terms will be studied. The typical magnitude of this L z a z term is estimated 

as 

^~^(-M 2 ~HT 8 eV. (A46) 
my so mif \XsoJ 

For example, the second order perturbation term, whose expression is similar to that in Eq. 
(1A39j) . contributes an amount of the order of 10~ 13 eV. In the higher order perturbation 
series the magnitude becomes even smaller and the contribution of the L z a z term to the 
S-T ± anticrossing gap would be of the order of 10~ 13 eV. 

Now we discuss the residual terms of the renormalized SO interaction in Eq. (IA41I) . 
These terms for two electrons can be rewritten in terms of the center-of-mass coordinates 
and the relative coordinates and again only the part associated with the relative coordinates 
contributes to the matrix element between the singlet state \S) and the triplet states |T±). 
The relevant part is given by 

H' - k x 

ren \ \ 

g [2{R^R n + ifr)pe + ±(J2^ + R 4V )U^ - 2(R 2 + ~£ 2 )p„ - R^U V - ^ 
v - ^2(R^R V + -^ri) Vv + l -(R^ + R^)U V - 2(R 2 V + ^ 2 )p 5 - R v r,U^ - ^ 



A 



The typical magnitude of this term is estimated as 

hH t h 2 f£ t \ 3 



(A47) 



mA| ml 2 V Ac 



~ 10~ n eV . (A48) 



In the higher order perturbation series, the contribution becomes much smaller. Thus the 
contribution from the renormalized SO interaction to the S-T± anticrossing gap is of the 
order of 10 -11 eV. 

Summarizing the above arguments within the linear-in-momentum SO coupling, we can 
conclude that the S-T-t anticrossing gap is mainly contributed by the SO induced Zeeman 
interaction H§° and this contribution can be eliminated in the first order by tuning the 
direction of the magnetic field, namely, by the magic angle tuning. However, the actual 
anticrossing gap is determined by other terms and the higher order perturbation terms of 
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H§° and is of the order of KT 12 ~ KT 11 eV. This magnitude is comparable to the hf 
interaction energy in the case of two delocalized electrons in a QD as discussed in Sec. IV. 

In general, the cubic-in-momentum spin-orbit term is present and more detailed argu- 
ments are necessary. The Dresselhaus SO term is originally given as 

V D = 7 A • a (A49) 
with A x = k x (k 2 y - k 2 z ) , A, = k y (k 2 z - k 2 x ) , A z = k z (k 2 x - k 2 y ) , (A50) 

where hki = p,i (i = x, y, z) including the vector potential due to a magnetic field. In the 
two-dimensional limit, we usually take the matrix element 

v£> = (<j>(z)\V D \<j>(z)) = ^{<j>(z)\k 2 M^)(-K<J x + k y a y ) (A51) 

with the ground state orbital <f)(z) in the ^-direction and put as 

= a D (p x a x - pydy) with a D = --(<j)(z)\k 2 z \(j)(z)) (A52) 

and call this the linear-in-momentum SO term which is already included in Eq. (J7J). Thus, 
in general, we can put 

V D = + V§ ] (A53) 

with V ( * ] = j{k x k 2 y a x - kyk 2 x o y ) , (A54) 

where vj^ 1 is called the cubic-in-momentum Dresselhaus SO coupling term. In the (£, 77, z) 
coordinate system this term is rewritten as 

V n = ^ UvlVn) - 2p>« + («,P€> - 2p?K] . (A55) 
After the unitary transformation in Eq. (1A26I) we have 

7 



2i 2i 

- r {-h{pe,p r ,} + i{{pt,p TI } - 2p 3 )£) + — (-h{pz, Pri } + i({p%pt:} - 2p\)rj) 



a. 



Ah 3 

= V§ ] + V§ )ren + • • • , (A56) 
where the results are given up to the first order in the expansion of Eq. (1A27|) and V^ ren is 

(3) (3) 

defined by the terms on the right hand side other than Vfj . Unfortunately, cannot be 
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eliminated. But the renormalized terms V^ ren are smaller in magnitude than the original 
terms due to the smallness parameter e = £t/^so ~ 10 -3 . Furthermore, those terms do not 
contribute to the matrix element between the singlet state \S) and the triplet states \T±) 
because they do not change the magnetic quantum number. Thus we have to consider the 

(3) 

effect of the original on the S-T± anticrossing gap. 

(3) 

The coupling for two electrons in the original coordinate system can be rewritten 
like in Eq. (|A9]) as 



y(V _ 7. 



{l U M + 1 U D +PxPyU y }E x + { Px (p 2 y + V) + -U x U yPy }a a 



1_ 



{\MpI + 4^) + p xPy U x }E y - {p y (pl + ^U 2 X ) + -U x U yPx }a y 



(A57) 



The operators T, x and S y do not change the magnitude of the total spin and thus they do 
not contribute to the matrix element (S\V^ \T±). Accordingly, we find 

(S\V£ } \T±) = TV2j^(g cm g rel \p x (p 2 y + -IlJ) + -U x U y p y \g cm e re i) 

(A58) 



1. 
4 



-i\f2 — (g cm g re i\py{p 2 x + -U 2 X ) + ^x^-yPx\9cme r el) 



1 
■r 



The linear terms with respect to the p x and p y can be included in Eq. (Al) to renormalize 
the linear-in-momentum SO interaction and the unitary transformation in Eq. (A26) can 
be redefined. Consequently, we have 



(S\V^\T ± ) = TV2jp(g rel \p x p 2 y =p % p y p 2 x \e rel ) . 
In order to estimate these matrix elements, we have to symmetrize the operators as 



1 



PxP v 



^{Px,Py} = ^(PxPy+PyPx) , PyP x ► ^{P V lPx} = 7>(P V Px + P xP V ) 



(A59) 



(A60) 



because p x and p y do not commute. The odd-parity excited orbital state of the relative 
motion can be approximated as 



1 



e rei( x ^v) = T^ x±i y) 9rei(x,y) with g re i{x,y) = - 
where the ^-coordinate part is omitted. Then we have 



r ^~ / 2 2 m 

■ e M-7^2( x +y)\ 



(9rel I ^ (PxP 2 y + p 2 y Px )\4el) = 1 ^ 1 



2£ 2 



1 



(9rel\^(PyP 2 x + p 2 x Py)\e^ el ) = ^~ 



h 3 



±^7T ± 



32£% 



+ 



=F 



\2U% 



2£ f V 2£ 2 32£% 8£ 2 



B 



B 



with 




(A61) 

(A62) 
(A63) 
(A64) 
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where £b is the magnetic length. The matrix element is calculated as 

(S\Vd 3) \ T +) = -^-^(9rel\^({Px,p 2 y} ~ i {Py,P 2 x })\e r el) 



for e ret (x,y) = e rel (x,y) 



(A65) 



In the same way, we find that (S\Vjp\TJ) = for e re i(x,y) = e^ el (x,y). In general, the 
S-T + and S-T_ crossing occur at different magnitudes of the magnetic field and we have to 
select an appropriate crossing point to realize the efficient electron-nuclear spin QST. In any 
case, by choosing the right S-T crossing point, we can eliminate the first order term of the 
cubic-in-momentum Dresselhaus SO coupling. 

(3) 

Now we consider the higher order perturbation terms with respect to Vp . The next 
third-order terms are given by 

, , v ^i^> 3) m)(ni^ 3) i^)(^i^ 3) i^) 

W 2^ (FATA - F(R'\\(F(TA - FATA ' 1 } 



S',T± 



(E(T±) - E(S'))(E(T ± ) - E{T' ± 

, (h] v (g|^ 3) |T4)(T4|^ 3) |To)(r |^ 3) |r ± ) 

^ 1 ] (E(T ± ) - E(T ))(E(T ± ) - E(T±)) 



(A67) 



where |S") and \T±) are the singlet and triplet states different from \S) and \T±). In the 
matrix elements of (T±\Vf) ,\Tq) and (Tq\V^\T±), the terms proportional to H x and H y in 
Eq. f]A57|) contribute. The general eigenstates can be written as 

|T+(T_)> = Vcm(R)0re*(r) 0(6) , (A68) 

|T ) = ^(R)0rd(r)-^(a(a)/3(6) + , (A69) 

where (p re i is an odd-parity excited state of the relative motion and ib cm and ^' cm are appro- 
priate orbital functions of the center-of-mass motion. The magnitude of the matrix elements 
in Eqs. (1A66I) and (1A67I) is typically of the order of 

(S|V®|T + )~£~l/ieV (A70) 



for 7 ~ 10 eV(A 3 )^ and £ t ~ 20 nm. Then the magnitude of terms in Eqs. (1A66|) and (1A67|) 
is estimated as 

1 /xeV x 1 fieV x 1 /xeV in _i2 



A A' 



~ 10~ 12 eV , (A71) 
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where A and A' are the orbital energy difference of the order of 1 meV. The higher order 

(3) 

perturbation terms with respect to become smaller furthermore. Finally, we discuss 
the contribution from Vp f n in Eq. f]A56|) which is smaller than by the smallness 
parameter e. The first order term (S\ V^ ren \T±) vanishes because V^ ren does not change 
the magnetic quantum number. In the higher order perturbation terms, V^ ren contributes 

(3) 

in combination with other interaction Hamiltonians, e.g., . For example, the second 
order perturbation terms have the magnitude of about 10~ 12 eV. The magnitude of higher 
order terms becomes smaller furthermore. 

Consequently, we can summarize that the S-T-t anticrossing gap due to the Dresselhaus 
and Rashba SO interactions would be of the order of 10~ 12 ~ 10~ n eV and is of the same 
order of magnitude as the hf coupling energy for the case of two delocalized electrons in a 
QD. This means that our protocol for the muclear spin quantum memory is not effective in 
the case of two delocalized electrons in a QD. 

So far we have considered the case of two delocalized electrons in a QD with an isotope 
atom of the host material. In the respect of the electron-nuclear spin QST time, the single 
electron charged QD with a donor impurity having the nuclear spin is more favorable. 
Especially favorable is the case where one electron is delocalized within the QD and the 
other electron is strongly bound to the ionized donor, as discussed in Sec. IV. In this 
case the hf coupling energy is of the order of 10~ 7 ~ 10 -6 eV. Now we examine the S- 
T± anticrossing gap due to the SO interactions. The most relevant term after the unitary 
transformation in Eq. ( 1A28I) is the SO induced Zeeman interaction H§°: 



H s z ° = g et i B 



t-ct£ + —a v \ B z - — B i + —B v ) a z 



(A72) 



This interaction for two electrons can be rewritten in terms of the total spin operators and 
the difference spin operators as given in Eq. ( 1A30I) and the matrix element between the 
singlet state \S) and the triplet state \T±) is given by Eq. (1A31|) . The orbital parts of the 
singlet and triplet states can be approximated by 

|vl> 5 ) = -L(0 a ( ri )0 b (r 2 ) + &(n)0 o (r 2 )) , (A73) 
|* T > = ^=[0 a ( ri )0 b (r 2 ) - <Mn)0 o (r 2 )] , (A74) 

where <p a (r) is the delocalized orbital in the QD and </>&(r) is the strongly localized donor- 
bound orbital. Hereafter the overlap integral between 4> a {r) and 0&(r) will be neglected 
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because the donor is doped in the peripheral region of the QD. Then we have 



(ttriei* s> = (*t|6 - &|*s> = (<f> a m a ) - <«I<M , (A75) 
(V T \v\Vs) = (^t\Vi ~ V^s) = {(paWa) - {<PbWb) ■ (A76) 

The center of the laterally symmetric confinement potential is chosen at the origin and the 
position of the donor atom is taken as (d, d, 0) in the (£, 77, z) coordinate system, leading to 

(*t|£i - 61*5) = -d , (V T \vi ~ m\*s) = -d . (A77) 

Then we obtain 

(T ± \H S Z °\S) = ±-- 



, ge^sBd 



s/2 



cos 02 _ smo9_ / smo9_ cos 09 . 

± l cos ^ 



(A78) 



\ A^ A^ 

This result suggests that the matrix element vanishes when 9 = tt/2 and tan<^_ = —X v /X^. 
However, this means the in-plane magnetic field and does not conform to the previous 
arguments developed for the case of a longitudinal (z-directed) magnetic field which is a 
prerequisite for the S-T level crossing. Thus the magic angle tuning is not possible here. 
Instead, in order to reduce the matrix element, we have to use the group IV elemental 
semiconductors in which the Dresselhaus SO terms are absent and the magnitude of the 
Rashba SO term might be reduced very much by careful tuning of the strain fields and 
the electric field. If we can achieve it/X^, P-t/\ ~ 10~ 4 , the above matrix element would 
be about 10~ 8 eV and be smaller than the hf coupling energy. At the same time, the 
contribution from the higher order perturbation series is much smaller. 

The arguments on the contribution from the renormalized SO interaction H ren in Eq. 
flA4ip can be developed in the same way as in the case of two delocalized electrons in a QD 
and the relevant matrix element is of the order of 10~ 12 ~ 1CT 11 eV. 

Finally, we examine the cubic-in-momentum Dresselhaus SO interaction Vp , although 
in the group IV elemental semiconductors this interaction is absent. As discussed before, 
Vj-, cannot be eliminated by the unitary transformation and should be considered in the 
original form. Vf) for two electrons is written as 



v {3) - -L- 



({Plx,Ply} + {p 2 x,P 2 y})^x - ({Ply,Pl x } + {p2y,P 2x })Zy 
+ ({Plx,Ply}-{P2x,ply})cTx-({Ply,pl x }-{p2y,pl x })cry , (A79) 
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where the momentum operators are symmetrized. Using the functions in Eqs. (I30I) and 
(]37p . we can show that 

(S\V£ ] \T±) = (A80) 

(3) 

by elementary calculations. The higher order perturbation terms with respect to Vp are 
present but the same arguments as in Eq. ( 1A71I) hold and their contribution is of the order 
of 10~ 12 eV. Finally, we discuss the contribution from V^ ren in Eq. flA56j) which is smaller 
than by the smallness parameter e. The same arguments below Eq. (1A71I) can be 
applied and the typical magnitude is estimated to be of the order of 10~ 12 eV. 

Consequently, in the case of the single electron charged QD doped with a donor impu- 
rity having the nuclear spin, the S-T± anticrossing gap due to the SO interactions can be 
neglected compared with the hf interaction energy. To achieve this situation, the group IV 
elemental semiconductors with the reduced Rashba SO interaction is favorable, which might 
be realized by careful tuning of the strain fields and the electric field. In these systems our 
protocol for the nuclear spin quantum memory would be fully effective. 
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APPENDIX B: INTEGRALS APPEARING IN THE EXPRESSION OF J 



The relevant integrals appearing in the expression of J in Eq. (1291) are given here. 

1 1 f°° 1 1 d 2 d 2 



1 /a/jVa 



2f, 2 ' v 2 2<r '- | 4( S 2 + | + 5 1 T ) 2£?" 

(B3) 



~/3V2 i 



=(0a(ri)06(r 2 )|-r-^ r |0 a (r 2 )0 fe (r 1 )) = i? a 



2 - 



2 2 

V, =<0«(r 1 )^ 5 (r 2 )|- j — ^ Mrjfofa)) ~ (0 a (n)| 6 ,.. |0 a (ri)) 

K|ri — r 2 | — ax | 

1 1 /-OO r2 

=R Wz -j=J Js (s 2 + l/ffi-^V + l/ffi^exp ^ + 1/£2) ] expha- 2 /^ 2 ] , 

(B5) 

_ la/3 1 / 2 , a 1 1 d 2 d 2 /£j-q 2 . . 

5= Vlp7 ( 2 + 2£f ) ( 2 + 2^ / eXp[ -W + 2(a + l/(^)) ] ' (B6) 

«= W 4, (B7) 



where R,aB,£t,&zi d,a and /3 are defined in Sec. Ill and the length variables £ t , i z and d are 
scaled by as to be dimensionless and the parameters q and s are also dimensionless. In the 
calculation, the transformation is employed for the Coulomb potential: 

1 1 f°° 

- = —= I ds exp[— s 2 r 2 ]. 



In the above, only the direct Coulomb term Vd is approximated by assuming that the QD 
confinement is much weaker than the donor atom confinement, i.e., it, iz ^ o-b/ \fo-i o,B/y/j3- 
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APPENDIX C: QUANTUM STATE TRANSFER BETWEEN AN ELECTRON 
AND A NUCLEAR SPIN IN THE GENERAL MIXED STATE 



Here we examine the quantum state transfer between an electron and a nuclear spin in 
the general mixed state, assuming 



Pn = 7l T)n n(T | + S\ |)„ „(| | + e\ |)„ n (| | + €* | |) n n (| 



(CI) 



where 7 and 5 are arbitrary non-negative real constants satisfying 7 + 5 = 1, e is an arbitrary 
complex constant and the suffix n attached to the bra and ket vectors indicates the nucleus. 
The electron spin is prepared in a pure state given by 



«l T)e + /3| l) t 



(C2) 



where a and j3 are arbitrary complex constants normalized as |a| 2 + |/3| 2 = 1 and the suffix e 
indicates the electron. Then the time evolution of the density matrix of the electron-nucleus 
coupled system is described by 



Pit) 



with p(0) = \ip) t 



\ a 2 7 


\a\ 


\a\ 2 e* 


\a\ 



a*(3j a*(3e |/3| 2 7 |/3| 2 e 
\a*(3e* a*{35 \/3\ 2 e* \{3\ 2 5 J 



where the bases of the matrix representation are arranged in the order of | |) e | f)n 
)n , I l) e | T)n an d | |) e | l) n . Then, using the expression of the hf interaction in Eq. 
have 

/ e -iAt/(4h) o o \ 

o e iAt/m cos m _ % e iAt/m sin m q 
-i e iAt/m sin At e iAt/m cos M 

e -iAt/(4h) 



(C3) 



(C4) 



Del I 

we 



e 



lit 



\ 



(C5) 



/ 
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and 



p(t = ttH/A) 



a(3*e \ 



( |a| 2 7 a/3*7 \a\ 2 

«73 7 |/3| 2 7 a*pe \(3\ 2 e 
\a\ 2 e* a(3*e* \a\ 2 5 a(3*5 



\ a* fie* 



12,* ~,* 



= (7l T)ee(T|+5| l)ee(l|+e| T>e e(l I + C*U>. e<T I) 
= Pe ® |^)n n(^| 
with |^) n = tt| |)n + /?| |)n , 



(C6) 
(C7) 



where p e is the same density matrix for the electron as in Eq. fICll) . This means that the 
QST or exchange of states between the electron and the nuclear spin is accomplished for the 
general mixed state of the nuclear spin. 



APPENDIX D: HYPERFINE COUPLING CONSTANT IN A LOCALIZED 
ELECTRON SYSTEM OF ZNSE:F 

The hf coupling between a donor electron and a donor nucleus is given by the Fermi 
contact interaction 2 ^ as discussed in Sec. IV: 

Vhf = ^9e^Bg n ^n S • I S(r - R) , (Dl) 

where S(I) is the dimensionless spin angular momentum operator for the donor electron 
(nucleus), /Ue(p; n ) the Bohr (nuclear) magneton, g e (g n ) the g-factor of the donor electron 
(nucleus), and r(R) denotes the position vector of the donor electron (nucleus). The donor- 
bound electron state can be represented by 

#(r) = F{v)u{r) (D2) 
with — / dr|u(r)| 2 = l, / dr\F(r)\ 2 = 1 , (D3) 

^0 Jvo J 

where u(r) is the Bloch function of the relevant conduction band normalized in the volume 
vq of a unit cell and F(r) is the envelope function. Then the hf coupling Hamiltonian is 
given by 

(Whf\*) = AS-I with A = ^-g e fi B g nf i n \F(R)u(R)\ 2 . (D4) 
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Now we infer the coupling constant A for the ZnSe:F in comparison with the case of Si:P 
based on Eq. (1D4|) . The squared amplitude of the envelope function can be roughly esti- 
mated by 

|F(R)| 2 ^ (D5) 
a B 

in terms of the donor electron Bohr radius given by 



= -7-5 , P6) 



m*e 2 



where m* is the isotropic effective mass of the donor electron and k is the dielectric constant 
of the host material. In the case of the anisotropic effective mass like in Si, Eq. (ID5I) is 
modified as 

|F(R)| 2 ^-1-, (D7) 

a t ai 

where a t m is the donor electron Bohr radius in the transverse (longitudinal) direction. Fur- 
thermore, it will be simply assumed that the Bloch function amplitude |w(R)| is not much 
different between Si and ZnSe. Then we can compare the coupling constant A between 
ZnSe:F and Si:P: 

^(ZnSe : F) = e (ZnSe)g n (F) a 2 (Si) a, (Si) 
A(Si:P) 9e(Si)g n (P) a|(ZnSe) ' 1 ] 

For Si, the g-factor of the donor electron is g e (Si) ~ 2.0, the nuclear g-factor of the P atom 
is <? n (P) = 1.13, and the donor electron Bohr radius is ~ 3.37(0.695) nm in the transverse 
(longitudinal) direction. On the other hand, for ZnSe, the nuclear g-factor of the F atom 
is (7 n (F) = 2.63, and the isotropic donor electron Bohr radius a# is ~ 3.01nm. The g- 
factor of the donor electron in ZnSe:F is not well known but it may be reasonable to guess 
that g e (ZnSe) ~ 2.0, because the measured g-factors of the donor-electron of many other 
impurities in ZnSe are accumulating around 2.0. Then the ratio of A is about 0.672 and we 
infer A(ZnSe : F) ~ 78.6 MHz from the value A(Si : P) ~ 117 MHz. 
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